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Abstract. The diffusion process through an ideal polymer network is studied by 
applying the effective medium theory (EMT) to a lattice-gas model. Polymers are 
modeled by random walks on the lattice bonds, across which molecules can hop with 
a certain probability. The steady state current of the system is calculated using the 
EMT and the results are compared to the simulations. 
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1. Introduction 

In this paper I apply the effective medium theory (EMT) to a prototype lattice- 
gas model simulating the permeation process in polymer networks. The permeation 
through polymer networks has been studied in a number of experiments and simulations. 
Experiments are usually performed on a thin polymer film subject to a constant gas 
pressure. The current is measured at the steady state and its relation to the temperature, 
the cooling rate, the size of molecules and so on are studied. [I] A prototype lattice-gas 
model, using random walks on lattice bonds to model the polymers, has been constructed 
in simulation studies. [2] This simple model is equivalent to a resistor network, in which 
the probability distribution of resistor on a bond is highly correlated to its neighbors. 
The correlations differ this model from the well-known randomly distributed resistor 
network and make it much harder to study. 

The EMT is a powerful tool to calculate the effective properties of random media. [3] 
It has been successfully used in calculating effective conductance of the completely 
randomly distributed resistor network when far away from percolation. In this case, 
the Bruggeman formula of the EMT provides sufficiently accurate results, |4] meaning 
that, when applied to our polymer systems, the Bruggeman formula yields excellent 
approximation when the polymer length is 1. When the polymer length gets longer, 
however, the error of the Bruggeman formula grows due to increasing correlations, and 
the theory needs be modified to account for the correlations. In this paper I develop 
a method to calculate the effective diffusivity, which gives fairly good results when the 
polymer length is relatively short (but greater than 1). Unfortunately, this method 
suffers from its fast growing complexity when the correlation increases, and thus is 
difficult to apply to very long polymer networks, or system close to percolation. 

2. Model and Theory 

Consider a d-dimensional lattice. We use coordinates x = (xi, Xd), where Xi, ■ ■ ■ ,Xd 
are integers, to label the lattice cells, and e^ {1 < i < d) to denote the base vectors. 
Polymers are randomly placed in the lattice, whose segments occupy the lattice bonds. 
In a numerical simulation the polymers are generated by random walks. A particle 
can occupy a lattice cell and can jump to one of its neighboring cells, across the 
polymer segments sitting between the two cells. The jumping probability depends on 
the number of polymer segments. Similar to the Ohm's law in a resistor system, the 
steady state current Ji(x), which flows from x + e^ to x, is related to the concentration 
drop Ei(x.) = c(x -|- ej) — c(x) by 



where a"jj(x) is the diffusivity. The diffusivity is a diagonal matric, whose diagonal 
elements o"m(x) are the jumping probability of particles across the bond separating cell 
X + ej and cell x. Throughout this section we adopt the Einstein summation convention 




(1) 
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that sums over repeated indices from 1 to d. The conservation of the number of molecules 
is written as 

Vi[a,,(x)A,c(x)] = 0, (2) 

where Aj is the forward finite difference such that, for any function /(x), Ai/(x) = /(x+ 
ej) — /(x), and Vj is the backward finite difference defined by Vi/(x) = /(x) — /(x — e,;). 

In a typical permeation experiment, a constant gradient of concentration Eq is 
applied on the system. It is convenient to define the effective diffusivity Sjj by the 
proportionality between the average current and the applied field. 

(J,) = S,, {Eo)y (3) 

where (■) denotes the spatial average. Assuming the system is isotropic, then we should 
have Sjj = a^Q^ij, where 6ij is the Kronecker delta. 

Let riij{x.) = (Jiji^x.) — T,ij. Eq. [2] can be solved using the Green function technique. 

c(x) = {Eo).x, + Y: G(x, x')V: (r^.,(x')A;.c(x')) , (4) 

x' 

where (^(x, x') is the Green function satisfying 

agg-V,A,G(x,x') = -5(x-x') (5) 
and vanishing on the boundary. Then E(x) is given by 

E,(x) = (Eo), + Yl G,,{^, x')r/,,(x')^fc(x') (6) 

x' 

where the d x d matrix G'(x, x') is given by 

G,,(x,x') = -A,A;.G(x,x'). (7) 

We also define susceptibility x(x) as a matrix such that 

Xij (x) {Eo)j = r]ij {^)Ej (x) , (8) 

which describe the field r7(x)E(x) as a response to the applied field Eq. Since 
(?7E) = (J - SE) = E ((Eo) - (E)) = 0, the average of susceptibility is zero, 

(X(x)) = 0. (9) 

Rewriting Eq. [6] in terms of x(x), we get the following Dyson equation 

X(x) = ^(x) + r^(x) Y ^(x, x')x(x'), (10) 

x' 

Decoposing the diagonal term 

X(x) = r/(x) + r/(x)G(x, x)x(x) + r]{^) ^ ^(x, x')x(x'). (11) 

and introducing 

X^(x) = [l-r^(x)G(x,x)]-Vx), (12) 
Eq. [11] can be simplified 

X(x) = x^(x) + x^(x) E G'(x, x')x(x'). (13) 

x^x' 
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Expanding the right hand side of the above equation and taking the average, we get 
(X(x)) = (x^(x))+ E (x^(x)G(x,xOx^(x') 

+ E E (x^(x)G(x,x')x^(x')G(x',x")x^(x")) + --- (14) 

X5^x' xVx" 

We classify the cells by the number of polymer segments sitting on its bonds. For 
example, when the dimension c? = 2, if each bond can only be occupied by one polymer 
segment at most, and if the probability of jumping across a polymer segment is q, then 
there are 4 different types of cells whose diffusivity matrices are 

...fi ?V3^fi (15) 



01/ \o q J \0 q 

respectively. Noticing that X'^(x) is determined by the type of cell at x, and in a 
homogeneous system, (^(x, x') only depends on x — x', Eq. [TUcan thus be rewritten in 
terms of statistics of different type of cells 



(X) = E /o(^) Xf + E E - ^-i) Xf ^(x' - x) xj 



i,j x't^x 



+ EE E /2(x'-x,x"-xV,J,fc)xfG(x'-x)xfG(x"-x')xf 

i,j,k x'^x x'^x" 

+■■■. (16) 

where /o(i) is the fraction of the type i cells, /i(x' — is the joint probability 

distribution of seeing a type i cell and a type j cell resting at x and x', respectively, and 
/2, /s, ■ ■ ■ are defined similarly. The function G'(x) = (5(0, x) can be easily calculated 
using the continuous Green function G'c(x, x') of the entire ci-dimensional space 

G,,(x) = A,V,G(0,x) (17) 



-1+^ , r-^H , .9^(0, xQ ^ / , / 1 



A,; 

Ixi- 

where 5{x) is the delta function. 

In principle, the probability distributions fi can be determined by doing statistics 
for a numerically generated polymer system. The computer resource required to do the 
statistics, however, increase exponentially with i and |x — x'|, |x' — x"|, etc. Therefore in 
practice we are only able to calculate the right hand side of Eq. [16] up to 3 terms, with 
the constrictions that |x — x'| < i?, |x' — x"| < i?, ■ ■ -, where > 1 is some constant 
integer we pick. Neglecting of the correlations outside the range R and higher order 
terms are justified by the fact that G'(x) decreases as jxj^'^. 

Combining Eq. [16] and Eq. [9] we get an equation that can be used to solve for the 
effective diffusivity Sjj. If only the first order term of the right hand side of Eq. [16] is 
kept, then it is only needed to solve (x^) = 0, which yields the classical EMT solution 
for the effective diffusivity. Noticing that G'jj(O) = —{l/d)5ij, the equation produces the 
well-known Bruggeman formula 

(l-p)(l-creff) ^ p{q - deff) ^ ^^^^ 



1 + {d- l)o-gff q + {d- l)cTgg 
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Table 1. Current measured in simulations Jgjjj^ and current Ji calculated by Eq. 
I16[ which includes up to i higher order terms, as a function of the polymer length /. 
Also shown in the brackets are the errors of each calculated current relative to the 
measurement. Note that Jq gives the classical EMT solution, which is independent of 
the polymer length. The probability to jump across an occupied bond is g = 0.1, and 
the occupation is p — 0.25. 



/ 




^0 


Ji 


J2 


^3 


1 


156.4 


155.7 (-0.45%) 


155.7 (-0.45%) 


155.9 (-0.32%) 


155.9 (-0.32%) 


2 


149.6 


155.7 (4.1%) 


151.5 (1.3%) 


151.3 (1.1%) 


151.1 (1.0%) 


4 


144.5 


155.7 (7.8%) 


150.5 (4.2%) 


148.9 (3.0%) 


148.7 (2.9%) 


8 


141.0 


155.7 (10.4%) 


149.7 (6.2%) 


147.1 (4.3%) 


146.7 (4.0%) 


16 


138.9 


155.7 (12.1%) 


149.4 (7.6%) 


145.8 (5.0%) 


145.4 (4.7%) 


32 


137.7 


155.7 (13.1%) 


150.1 (9.0%) 


146.0 (6.0%) 


145.5 (5.7%) 


64 


137.2 


155.7 (13.5%) 


150.6 (9.8%) 


146.1 (6.5%) 


145.6 (6.1%) 



where p is the fraction of occupied bonds in the system. 

In a previous study it has been shown that for a system with the polymer length 
i = 1, Eq. [T8]gives an excellent approximation to crgg.|3] However, if the polymer length 
i > 1, the higher order terms, which describe the correlation between cells, significantly 
revise the results. It is also worthy of noting that this method of calculating the effective 
diffusivity works best when the system is far away from the percolation threshold, which 
means either the temperature has to be high, or the mass density be very low. Only in 
this regime the correlation is not crucial, and this method can be a good approximation. 

For long polymers or system not so far away from percolation, or, in other word, 
when correlations are important, one could improve the result of the EMT by including 
more higher order terms in Eq. [161 This is illustrated by our calculation result shown 
in Tabled], where we list the current at g = 0.1 and p = 0.25 calculated through Eq. 
[T6l using up to 3 higher order terms in the expansion, in comparison with both the 
current measured by numerical simulations (for details of the simulation, see [5]) and 
the current calculated by the classical EMT approximation. The correlation range R 
used in the calculation is 3. Only single occupation on bonds is allow in the calculation, 
because multiple occupation produces too many types of cells to do statistics for i? > 1 
or include more than one higher order terms. 

The results clearly show that the classical EMT yields extremely accurate result 
for i = 1, but produces considerable errors when £ > 1. Taking account of correlations 
by including the higher order terms can significantly improve the calculation. 

3. Summary and discussion 

In this paper I have used the EMT to calculate the effective diffusivity of a lattice-gas 
model simulating permeation through polymer networks. Such a model is equivalent to 
resistor network, which is described by a Laplace Equation at the steady state. The 
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solution of this equation can be expressed by a Dyson expansion. One needs to include 
higher order terms in the expansion to account for the correlations of the distribution 
of polymer segments. When system is far away from percolation, we have calculate 
the steady state current using up to 4 terms in the expansion and get reasonably good 
result compared to simulations. The error of the calculation increases with the polymer 
length, since longer polymers have stronger correlations. 

From Table [1] one can see that the current is decreasing function of the polymer 
length. It can be qualitatively explained by considering the blob (the set of cells 
that carry current at g = 0) size, which also decrease with the polymer length. [5] A 
possible future work would be using the EMT to provide some quantitative insight to 
this phenomenon. 
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